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ABSTRACT 

We analyze the clustering properties of ultraviolet selected galaxies by using GALEX - SDSS 
data at z < 0.6 and CFHTLS deep u' imaging at z ~ 1. These datasets provide a unique 
basis at z < 1 which can be directly compared with high redshift samples built with similar 
selection criteria. We discuss the dependence of the correlation function parameters (r , S) on the 
ultraviolet luminosity as well as the linear bias evolution. We find that the bias parameter shows 
a gradual decline from high (b s > 2) to low redshift (b s ~ 0.791"'^). When accounting for the 
fraction of the star formation activity enclosed in the different samples, our results suggest that 
the bulk of star formation migrated from high mass dark matter halos at z > 2 (10 12 < M m i n < 
10 13 M Q , located in high density regions), to less massive halos at low redshift (M min < 10 12 M Q , 
located in low density regions). This result extends the "downsizing" picture (shift of the star 
formation activity from high stellar mass systems at high z to low stellar mass at low z) to the 
dark matter distribution. 



Subject headings: Galaxies: UV - Correlation Function Evolution - Star Formation - Downsizing 
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1. Introduction 

Accumulated evidence shows that the cosmic 
Star Formation Rate (SFR) has been decreas- 
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in g from z ~ 1 by a dramatic factor of about 
5 iHopkinsI l2004t iLillv et all Il996t [Madau et al ' 



Schiminovich et al.1 



Wilson et all 200^ 



2005; Sull ivan et al 
This is linked to the 
decrease of the contribution of the faint galaxies 
that dominate the star formation density, and 
to the strong decline of the most ultraviolet- 
luminous galaxies with time, given the rcdshift 
evolution of the 15 00 A luminosity function 
(jArnouts et al. 2005t) . Another as pect of this 
evolu tion, known as "downsizing" (jCowie et al 
is the observation that star formation ac- 
tivity shifts wi th time from high to low stellar 
mass systems (iBundv et al.l 2005 : Jimenez et al 



20051: 1 Juneau et alJl2005tlHeavens et al.ll2004 and 
references therein). 

The star formation history results from the in- 
terplay between the physical processes driving the 
star formation fueling (gas cooling) and regula- 
tion (feedback), both closely related to galaxy en- 
vironment. Recent simulations show that about 
half of the galaxy gas is accreted through a cold 
mode, which dominates at high redshift in high 
density environments, and shifts to low density 



environments in the local Universe (Keres et al 



2005) . The type of the dominant feedback process 



is expected to depend on galaxy host halo ma ss: 
supernovae explosions (e.g. Benson et alJ 2003 ) at 
low mass, an d AGN (e.g. ICroton et al.l 1200(1 ) at 



high mass. iCattaneo et al.l (|2006h show that the 



IS 



introduction of a critical halo mass above which 
there is a complete shutdown of cooling and star 
formation is efficient to reproduce the bimodality 
in g alaxy properties obs erved in the local Universe 
(e.g. iBaldrv et al"]l2004h . 

In this paper, we propose to set constraints 
on the roles of these different processes through 
cosmic time by assessing the spatial distribu- 
tion of star formation in the Universe from 
high to low redshifts. A convenient method 
is to study the clustering properties of rest- 
frame ultraviolet (UV) selected galaxies. This 
has already been performed at high redshifts 
using Lyman Break Galaxies (LBGs) samples 
to show that, at these epochs, star formation 
is highly clu stered and concentrated in over- 
dense regions ( Adelberger et al.l 20051: Allen et al 



2005: Arnouts et all l2002t iFoucaud et all 



20031: 



Giavalisco k. Dickinson 200ll ). The study of the 



actively star forming galaxies has now been made 
possible in a homogeneous way with the combina- 
tion of restframe UV data collected from z ~ 4 to 
z = 0. To extend high-z studies, we use GALEX 
observations in the recent Universe and CFHTLS 
deep imaging at z = 1. Wc compute the angular 
correlation function (ACF) of star forming galax- 
ies and derive their bias and its evolution. 



In a companion paper, iMilliard et alJ (2007, 
hereafter Paper I), we describe in detail the 
methodology and the first results of the angu- 
lar correlation function measurements of UV se- 
lected galaxies using a GALEX sample at z < 0.6. 
Section [2] summarizes the sample properties and 
presents a new restframe UV-selected sample from 
the v! band deep CFHTLS imaging survey that 
we use to extend the analysis to higher redshift 
(z ~ 1). We then investigate the dependence on 
redshift and UV luminosity of the clustering prop- 
erties: ro, S in sect. 02 bias in sect. HI In the last 
section we discuss the evolution of the preferred 
sites of star formation over the last 90% of the age 
of the Universe. 

All magnitudes have been corrected for Galac- 
tic extinction using the E(B-V) v alue from the 
ISchlegel et al. ( 19981 ) maps and the ICardelli et al. 
(jl989l ) extinction law. Throughout the paper, 



redshift evolution of the clustering properties of 



we adopt the following cosmological parameters: 
n m = 0.3, fl A = 0.7, H = 70 km s^Mpc -1 . 

2. Samples description 
2.1. GALEX 

In this work, we use the same subsample of 
GALEX Release 2 (GR2) Medium Imaging Sur- 
vey (MIS) fields cross-matched with SDSS DR5 
presented in Paper I, and we refer to this paper 
for a full description. We recall here the main 
characteristics of the selection. We only keep 
GALEX objects with SDSS counterparts within 
a search radius of 4" and use the closest SDSS 
match. We select galaxies as objects with SDSS 
type equal to 3. We use the half of the MIS fields 
from the GR2 dataset with the lowest Galactic ex- 
tinction {{E(B - V)) < 0.04). Photometric red- 
shifts are computed using an empirical method 
()Connollv et al.lll995h trained on SDSS spectro- 
scopic counterparts. The standard deviation es- 
timated from the SDSS spectroscopic redshifts is 
a = 0.03. We then use a template fitting proce- 
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Fig. 1. — NUV absolute magnitude-photometric red- 
shift relation in the GALEX sample. The colors code 
the type according to the best fitting template: red 
represent elliptical types, green spirals, and blue irreg- 
ulars. The dashed lines indicate the additional cuts 
adopted: z phot < 0.6 and -21.5 < NUV ahs < -14. 
The same cuts hold for FUV. 



dure (jArnouts fc Ilbert. Ilin preparation! ) to derive 
UV luminosities. Our starting samples include ob- 
jects with FUV < 22 or NUV < 22. 

The NUV absolute magnitude vs photomet- 
ric redshift relation is shown in figure [TJ The 
colors code the galaxy type determined using a 
SED template fitting procedure: red represent 
elliptical types, green spiral and blue irregular. 
We restricted hereafter the samples to —21.5 < 
NUV^s < -14. and 0. < z phot < 0.6 (dashed 
lines on fig[T]). The same cuts have been applied 
to the FUV sample. 

In the following, we consider both FUV and 
NUV bands and we divide the samples in two bins 
according to the mean UV absolute magnitude. 
The figure [2] shows the photometric redshift distri- 
butions of the GALEX samples; the table [T] sum- 
marizes the properties of the samples. 

2.2. CFHTLS 

The CFHTLS-Deep survey consists of deep 
multi-colour images collected through the u'g'r'i'z' 
filters over four independent areas of 1 deg 2 each 



Fig. 2. — Redshift distribution of the subsamples cut 
in absolute UV magnitude: M < (M), solid lines; 
M > (M), dashed lines, FUV is shown as blue and 
NUV as red. 



and reaching the limiting magnitude of i' AB ~ 26. 
In this work, we use the official CFHTLS data 
release T0003. For a full pre sentation of the 



CFHTLS-Deep survey, we refer to lSchultheis et al 



( 2006l f1. We built specific masks from the u-band 
images to mask out stars, chips edges' and ar- 
tifacts. The total solid angle of the four fields 
used after masking is 3.1 deg 2 . The star/galaxy 
separation is ba s ed on the same method as 
McCracken et al.l (|2003h with the half-light ra- 
dius versus u magnitude plot. This selection has 
been applied down to u = 23. Beyond this limit, 
we combine the photometric criterion with the 
star/galaxy classification deriv ed from the photo 



metric redshift code Le Phare, (jArnouts fc Ilbert. 
in preparation! ) . 



To construct the sample of UV selected galaxies 
at z ~ 1, we adopt a u' magnitude limit of u' — 24, 
which ensures a genuine C/U-selected sample as 
the it' effective wavelength (3587 A) corresponds 
to 1848 A at our mean redshift (z) = 0.94. The 
fraction of objects lost (without any redshift selec- 



1 see also http://www.cfht.hawaii.edu/Science/CFHLS/ 
and http : //www. ast . obs-mip . f r/article204 .html 
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Fig. 3. — Redshift selection based on the (g' — r') vs 
fr'-i') diagram for the CFHTLS-D1 field. The small 
dots show galaxies with z p hot < 0.7 and big dots galax- 
ies with Zphot > 0.7. The line represents the adopted 
color-color selection criterion. 



tion) with a i' — 26 cut is on average 0.07% over 
the four fields at u' — 24. The redshift selection 
of the sample is based first on a color-color selec- 
tion and then on the photometric redshift selec- 
tion. We do not adopt a single selection based on 
the photometric redshifts because of the variable 
accuracy of the method due to inhomogeneous ex- 
posure times in the five bands for the different 
fields. 

First, we use a color-color selection to isolate 
galaxies with z > 0.7, based on VVDS photo- 
metri c redshifts estima tion, relying on multi-color 
data ( Ilbert et al.ll2006h . As shown in figure [H the 
(g — r) versus (r — i) selection criterion is efficient 
to separate galaxies at z > 0.7 (big dots) from the 
lower redshift population (small dots). The line 
shows our separation criterion. 96% of galaxies 
with Zphot > 0.7 are located below the line while 
less than 10% of low z objects (z phot < 0.7) fall in 
the same region. 



The photometric redshifts are computed by us- 
ing Le Phar e code and by ado pting the method 
described bv lllbert et al. ( 2006h . The comparison 
with the spectroscopic redshifts, obtained by the 



WPS in the best pho tometric field (CFHTLS-D1, 
Le Fevre et"aH (|2005t )). for our u 1 selected sample 
shows an accuracy of a(Az/(l + z)) = 0.03 with 
4% of outliers (defined as Az > 0.15 x (1 + z)). 

In figure [4] we show the photometric redshift 
distribution of the galaxies selected with the color 
criterion (dashed histogram). The final sample is 
obtained by further selecting objects with 0.7 < 
Zphot < 1-3 (solid histogram). The absolute mag- 
nitudes in the GALEX bands are derived from the 
best fitting SEDs whose NUV-rest flux are well 
constrained by the u', g' , and r' bands in the red- 
shifts range (0.7 < z < 1.3). Note that as the v! 
filter shifts to FUV wavelengths at z ~ 1, absolute 
magnitudes depend very weakly on k-correction 
and best-fit fitting SEDs. As for GALEX sam- 
ples, we divide the CFHTLS sample in two bins 
according to the mean FUV absolute magnitude 
and the resulting redshift distributions are shown 
in Fig[U 

The global properties of the CFHTLS UV samples 
are given in table [2] 

3. Redshift evolution of the correlation 
function of UV-selected galaxies 



We compute the ACF using the lLandv fc Szalav 
(|l993l ) estimator. We assume that the ACF is well 
approximated by a power-law: w(6) = A w 6~ s ; we 
use a variable Integral Constraint (IC) with S as 
free parameter during the power-law fitting pro- 
cess, an d estimate the I C wit h the same method 
used by Roche &; Ealesl ( 19991 ). We derive correla- 
tion lengt hs (rn) for eac h sample from the Limber 
equation (IPeeblesll^Soh . using the corresponding 
redshift distribution. These quantities, as well as 
the bias parameteiQ, are summarized in table [T] 
and table [2] The effects of the dust internal to 
galaxies have again been neglected. 

In figures [5] and [H we show the ACFs of the 
GALEX and CFHTLS samples respectively. The 
ACFs are derived for the global samples and for 
two sub-samples with UV absolute luminosity 
brighter and fainter than the mean < UV a b s > 
of each sample. The angular scales probed for the 
GALEX samples are 0.005° to 0.4° (correspond- 
ing respectively to comoving distances 0.07 Mpc 
and 5.7 Mpc at z = 0.2), while 0.002° to 0.4° 



See Paper I for details on the computations. 
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Fig. 4.— Redshift distribution of the CFHTLS 
sample: the dashed histogram shows the photomet- 
ric redshift distribution of the galaxies selected with 
the color criterion alone, while the solid histogram 
shows the final redshift distribution after selecting ob- 
jects with 0.7 < Zphot < 1-3. The dotted and dot- 
dashed histograms show the redshift distributions of 
the FUVabs < -19.41 and FUVabs > -19.41 respec- 
tively. 



for the CFHTLS samples (resp. 0.11 Mpc and 
23 Mpc at z = 1). These ACFs are fairly well 
fitted by power-laws, even if a small dip appears 
at small scales in the FUV GALEX samples and 
also in the CFHTLS bright one (see sec. I3~2~2"j) . 
The higher surface density of UV-selected galaxies 
at z ~ 1 allows a less noisy estimation of the ACF 
at these epochs than at z < 0.4. 

3.1. Clustering segregation with FUV lu- 
minosity 

The dependence of tq on FUV luminosity in 
GALEX and CFHTLS samples is shown in fig. 
along with results from higher redshift studies 
(z > 2^1 As the different surveys probe differ- 



3 We choose the FUV absolute magnitude for the compar- 
ison as most of high redshift samples are FUV restframe 
selected, and GALEX results are not strongly dependent 
of the UV band. 

The mean absolute magnitudes of the LBG samples have 
been obtained by deriving an average apparent magnitude 



Fig. 6. — Angular correlation function of the 
CFHTLS subsamples cut in absolute UV magnitude: 
M < (M), squares; M > (M), triangles and for com- 
parison the total sample (circles). The curves show 
the best fit power-laws with the Integral Constraint 
correction terms subtracted. The upper axis shows 
the comoving scales corresponding to angular scales 
at z = 0.9. 



ent parts of the UV luminosity function with lit- 
tle overlap, it is difficult to draw firm conclusions. 
Nevertheless, significant differences between the 
samples are apparent: 

• We use as reference the correlation func- 



tion results from j Adclbcrg er et al 
Allen et all (120051). Foucaud et al 



(2005), 



2003) 



Giavalisco fc Dickinson (1200 II) . iLee et al 
(|2006h and lOuchi et alj (|200lh at z > 2. At 
these redshifts, all studies conclude a signif- 
icant segregation of ro with UV luminosity 
(the more luminous the more clustered) in 
the range -23 < FUV^s < -20. 



At 



1, a positive correlation of rg with 



FUV is still observed for -20 < FW abs < 
— 19. Notably, our value of tq at FUV a hs ~ 



from the galax y counts, an d assuming a k-correction of 
2.51og(l + z). lOuchi et all l|2005l ) do not provide their 
counts, so we computed the expected mean absolute mag- 
nitude given their li miting absolute magnitude a nd the lu- 
minosity function of Sawicki & Thompson (2006). 
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Table 1 

galex samples description, power-law best fits parameters, comoving correlation lengths and 

BIAS. 





FUV samples 




NUV samples 


All 


FUV abs < -18.3 F!7V abs > -18.3 


All 


iVC/Vkbs < -18-8 WV aba > -18.8 



N gal * 


42065 


22082 


19983 


97038 


52567 


44471 


<Ft/y aba ) 


-18.3 


-18.96 


-17.57 


-18.23 


-18.76 


-17.61 


aFUV a.bs 


0.91 


0.52 


0.67 


0.96 


0.72 


0.82 


(ivt/v abs ) 


-18.58 


-19.15 


-17.95 


-18.8 


-19.43 


-18.05 


CjV(7V abs 


0.84 


0.52 


0.67 


0.91 


0.44 


0.69 


<*} + 


0.18 


0.23 


0.13 


0.25 


0.32 


0.17 




0.08 


0.07 


0.05 


0.11 


0.09 


0.06 


n gal [10- 2 Mpc- 3 ] 


2.78±1.03 


0.15±0.03 


2.54±0.95 


2.16±1.19 


0.14 ±0.03 


1.95± 1.06 


A w x 10 3 


9.4+- 


7 3+ 34 

' -°-2.2 


„„ q + 17.4 


3 6+°- 6 


4 l + 1 -° 


22 4+ 5 ' 7 


<5 


0.74±0.05 


0.79±0.1 


0.48±0.09 


0.86±0.04 


0.87±0.05 


0.52±0.05 


r [Mpc] 


4 6+ 06 


4 6+ 09 


5.4tl; 5 


4 1+0-3 
*' A -0.3 




c c+0.8 
°-°-0.7 


6s 


74+o.os 

u - '^—0.07 


0.76±g i 3 


0.83+g;S 


n 69+ 05 


0.83+°;°a 7 7 


n »k+ 09 



* Number of galaxies in the samples 
^according to photometric rcdshifts 



Note. — The amplitude and slope of best fit power laws to the angular correlation function, and hence the comoving correlation length 
account for the Integral Constraint correction, as described in Paper I. 



Table 2 

CFHTLS SAMPLES DESCRIPTION, POWER-LAW BEST FITS PARAMETERS, COMOVING CORRELATION LENGTHS AND 

BIAS 





All 


FUV^s < -19.41 


FUV^s > -19.41 


Ngal* 


17098 


8507 


8591 


(FUV ahB ) 


-19.41 


-19.89 


-18.94 


a FUV abB 


0.6 


0.34 


0.36 


(AW abs ) 


-19.81 


-20.19 


-19.43 


CT iWV aba 


0.53 


0.4 


0.33 


(z)t 


0.94 


1.04 


0.84 




0.16 


0.15 


0.09 


n gal [l(r 3 Mpc- 3 ] 


3.27±2.90 


0.62±0.41 


2.66±2.58 


A w x 10 3 


2 7+ 1 - 4 


9 H+ 1 - 2 


q q+1.0 

°-°-0.7 


8 


0.7 ±0.09 


0.74 ±0.09 


0.76 ±0.07 


ro [Mpc] 


4.92 +g;| 


5.48l° ;| 


4.66I 023 


68 


^ -0.07 


i no+0.06 
1 -°°-0.07 


1 16+ 003 

1 - ±D -0.04 



'Number of galaxies in the samples 
t according to photometric redshifts 

Note. — Same note as in table [T] 
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Fig. 5. — Angular correlation function of the GALEX subsamples cut in absolute UV magnitude: M < (M), squares; 
M > (M), triangles and for comparison the total sample (circles). Left panel, FUV; right panel, NUV. The curves 
show the best fit power-law not corrected for the Integral Constraint bias. The upper axis shows the comoving scales 
corresponding to angular scales at z = 0.18 (left) or z — 0.24 (right). 



—20 is very close to that of lAdelberger et al 
(|2005l ) obtained from z ~ 2 samples. 



• At z < 0.3, we probe a fainter luminosity 
range (—19 < FUV a bs < — 17), and a weak 
anti-correlation of ro with FUV is apparent, 
though given the error bars, it is compatible 
with no FUV luminosity segregation of ro. 

The values of ro as a function of FUV lumi- 
nosity for FUV selected samples at different red- 
shifts follow a unique smooth curve, with a sig- 
nificant slope at the bright end FUV^hs < —19 
and a flat or slightly negative slope at the faint 
end. A similar segregation is observed with B 



luminosity at low redshift, 


with optical selec- 


tion criteria (|Benoist et al. 


19961 Guzzo et al. 


2000; Norberget al. 20021 


Willmer et al.1 Il998t 


Zehavi et al. 2005). In particular. iNorberg et al.l 


(2001) showed that for blue-selected galaxies r n 



increases only slowly for galaxies fainter than 
Lf , while it varies strongly for galaxies brighter 
than . I ndeed , using TV-body simulations 
Benson et al . (2001) showed that Lf could be a 
natural boundary in the distribution of the halos 
hosting galaxies, galaxies fainter than Lf being 
hosted by a mix of low and high mass halos, while 
galaxies brighter than if hosted by more and 
more massive haloes. To check if FUV* could play 
a similar role in UV samples, we show in figure[8]ro 
as a function of (FUV a hs) ~ FUV*, where the evo- 



lution of FUV* with z is tak en from Arnouts et al 
(|2005[) (for z < 1) and ISawicki fc Thompsor 
(|2006l ) (for z > 1). The luminosity dependence of 
ro changes noticeably when expressed as a func- 
tion of FUV a hs — FUV*, as two different trends 
are observed according to the redshift range: 

• At z > 1, for the high z samples and our 



7 



redshift redshift 




10 



Adelberger et al (2005) □ 
Arnouts et al (2002) * 
Foucaud et al (2003) a 
Giavalisco & Dickinson (2001) o 
Heinis et al (2004) » 
GALEX FUV • 
GALEX NUV ■ 
CFHTLS u » 




-23 



-22 



-21 



-20 
< FUV ab 



-19 



-17 



10 




Adelberger et al (2005) □ 
Arnouts et al (2002) * 
Foucaud et al (2003) a 
Giavalisco & Dickinson (2001) o 
Heinis et al (2004) » 
GALEX FUV • 
GALEX NUV - 
CFHTLS u * 



<FUV abs )-FUV, 



Fig. 7. — Dependence on absolute FUV magnitude 
of the correlation length ro for low and high red- 
shift restframe {TV-sele cted galaxies: open squares, 
lAdelbereer et all lj2005l); open star , lArnouts et al] 
l|20(% open triangles. [ Foucaud et al] (|2003l ) ; open cir- 
cles, iGiavalisco fc Dickinson (2001). Our results are 
presented as filled circles (FUV) and filled squares 
(NUV) for the GALEX samples and as filled stars for 
the CFHTLS samples. The mean redshifts of the sam- 
ples are color-coded. Note that the results from a given 
study (including ours) are not all obtained from inde- 
pendent samples. Hence we distinguish global samples 
by plotting them with a bigger symbol size, in this fig- 
ure and in the following ones as well, unless otherwise 
stated. 



CFHTLS sample, the behavior of ro with 
FUV a bs — FUV* is qualitatively compatible 
with the monotonic trend described above, 
the brighter galaxies being more clustered. 

• At z < 0.5 (GALEX samples) a radically 
different behavior of ro vs FUV a \>s ~ FUV* 
is seen. An anti-correlation or no correla- 
tion (given the errorbars) is observed, with 
brighter samples showing slightly lower ro 
than fainter ones. 

This suggests that the luminosity segregation 
mechanisms of the clustering at low redshifts work 
in a different regime, or that FUV* is not the rel- 
evant variable. 



Fig. 8. 
FUV*. 



Same as in figure [7] but versus FUV a h B 



3.2. ACF slope segregation with FUV lu- 
minosity 

3.2.1. ACF slope 

The slope of the ACF (5), which describes 
the balance between small and large scale sep- 
arations, is an important indicator on the nature 
of the spatial distribution of a given population. 
In paper I, we found that the estimates of the 
slope inferred from the global GALEX samples 
(S ~ 0.81 ± 0.07) are steeper than those de- 
rived from optically selecte d blue galaxies in the 
local Universe: S ~ 0.6 (Budavari et alj l2003t 
Madgwick et alj|2002l: IZehavi et alj|2002l l2005h . 
In figure O we now analyse the dependence of 
S on UV luminosity for the different samples 
(GALEX samples: filled circles for FUV and filled 
squares for NUV; CFHTLS samples: filled trian- 
gles; high z samples: empty squares, open circles 
and crosses; all points color-coded with redshift). 

At z > 3, the compilation of measurements 
show ed here, and espec ially those at at z = 4, 



from Ouchietal 



( 20051 ) indicate that the ACF 
slope steepens at higher UV luminosities. 

( 20051) claimed that this trend, 



Ouchi et al 



well modeled in the Halo Occupation Distribu- 
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Fig. 9. — Dependence on absolute FUV magnitude 
of the slope 8 of the correlation function for low and 
high redshift restframe ?7V-selected galaxies. Results 
shown here come from stu dies allowing 5 a s a fre e 
parameter: open squares. | Adelberger et all ()2005l ): 
open circles , Giavalisco fc Dickinson! I 200lf) ; crosses. 
lOuchi et all {2005) (results from the fit with Integral 
Constraint). Our results are presented as filled circles 
(FUV) and filled squares (NUV) for the GALEX sam- 
ples and as filled stars for the CFHTLS samples; the 
error on S comes from the fitting procedure and the 
horizontal bars for our samples reflect the standard 
deviations of FUVa.be- The mean redshifts of the sam- 
ples are color-coded. 



tion (HOD) framework, is not an actual slope 
variation but it is related to the halo occupancy. 
Based on HOD models, they show that the con- 
tribu tion of satelli t e gala xies ('one halo term', see 
e.g. 



Zehavi et al 



(2004)) increases when select- 
ing brighter galaxies, by enhancing the small scale 
signal of the ACF (r < 0.35 Mpc). This effect 
produces an apparent steepening of the observed 
slope S. 

Our GALEX sample at z < 0.4 seems to pro- 
duce a similar although less pronounced effect. 
Our current low z GALEX data do not allow to 
perform detailed comparison between observations 
and HOD models, but we have investigated if the 
observed steepening with luminosity can be par- 
tially due to the small scale component. We fitted 
the GALEX ACF only at scales r > 0.4 Mpc (see 



sec. I3.2.2|) or r > 0.7 Mpc in order not to in- 
clude the one halo term component. We do not 
observe any significant change with respect to our 
initial slopes. However doing so we face at large 
scales the problems of lower signal-to-noise ratio 
and more important contribution of the Integral 
Constraint bias that prevent us to make firm state- 
ments. This test thus relies on the efficiency of our 
power law fitting process in recovering the true 
ACF (see Paper I). 

In other words, at low redshift, we do not see 
evidence that the one halo term plays a major 
role in the slope of the ACF, as observed at 
high redshift, which i s expected from simulations 



(jKravtsov et al.ll2004f ). Hence this indicates that 
our clustering parameters (rn, S and bias, b$) 
reflect the large scale clustering of star forming 
galaxies, which enables us to make comparisons 
with analytical predictions for the clustering of 
Dark Matter Haloes. 

3.2.2. Dip in the ACF ? 

The ACFs derived for the various GALEX and 
CFHTLS samples are globally well described by 
a power-law, but some of our ACFs show a lit- 
tle dip around 0.35 Mpc, the GALEX FUV ones 
for instance, and also the brightest CFHTLS sam- 
ple at z ~ 1 (at a slightly larger scale ~ 0.5 
Mpc). This recalls the departure to the power- 
law observed in other surveys and interpreted 
as the transition between the one and two halo 



terms in the HOD framework. IZehavi et al.l ([2004) 
showed that this transition occurs at ~ 1.5 — 3 
Mpc for r-band selected galaxies. It is expected 
that this scale should be shorter for bluer galax- 
ies, residing in less mass i ve hal os, as showed by 
Magliocchetti &: Porcianil (|2003l ) in observations 



and lBerlind et al.l (|2003l . see their fig. 22) in simu- 
lations, with a transition scale for late-types galax- 
ies at ~ 0.45 Mpc, close t o what we observe. Fi- 
nally and very interestinglv lOuchi et al. | (Hooi) ob- 
serve this transition for LBGs at z ~ 4 at 0.35 
Mpc, the same comoving scale as we get. 

Comparing measurements with predictions 
from HOD models is a natural perspective of this 
work, to probe the redshift evolution of the halo 
occupancy of star-forming galaxies. This will be 
addressed in details in a forthcoming paper with 
enlarged datasets. 
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4. Bias of star-forming galaxies from z = 4 

to z = 

The link between the properties of the galaxy 
distribution and the underlying Dark Matter den- 
sity field can be accessed via the bias formal- 
ism. The bias parameter is indicative of the 
masses of the dark matter halos that preferen- 
tially host the observed galaxy population (e.g . 
Giavalisco fc Dickinson! |200lt IMo fc Whitel 120021: 
Ouchi et al J 120041 ). i.e. in our case, actively star 
forming galaxies. The DM h alo bias is a direct 
output of IMo fc White] (|2002f ) models. For galax- 
ies, we assume a linear bias to convert ro in <7s.g 5 
a c ommon, though quest ionable, assumption (see 
e.g. iMarinoni et al.ll2005T ). 



pies discussed above The bias values for our 
GALEX and CFHTLS samples are reported in 
tables Q] and [2] respectively. 

The observed bias of star-forming galaxies shows 

a gradual increase with look ba ck time: at z > 2, 

UV galaxies are strongly bi ased (jGiavalisco fc Dickinson 
20011 iFoucaud et alll2003h . with b 8 > 2, and at a 



given redshift the bias increases with FUV lumi- 
nosity (FUV luminosity segregation). At z ~ 1, 
the mean bias is {b 8 ) = 1.26 ± 0.06, indicating 
that star-forming galaxies are closer tracers of 
the underlying mass distribution at that time. 
At z < 0.4, given the error bars, the mean bias 
is consistent with 0.8 for all GALEX samples 
({bs) = 0.79lo;J 8 ), a slight anti-bias independent 
of the UV luminosity. 
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Fig. 10. — Evolution with redshift of the bias of rest- 
frame [TV-selected galaxies (symbols as figure[7]). The 
curves show the effective bias of halos more massive 
than Mmin (color-coded) according to IMo fc Whitel 
(2002) (see also Paper I). The inset is an enlargement 
of the low redshift area. 



In figure [TOl we also s how the effective bias evo- 
lution derived from the IMo fc Whitel ([2002) for- 
malism for different minimum Dark Matter Halo 
(DMH) mass thresholds. A comparison can be 
made to the bias of star forming galaxies, if one 
assumes that most haloes do not host more than 
one star-forming galaxy. This coarse assumption is 
likely inaccurate for star-forming galaxies selected 
at high redshif ts in FUV with a well developed 
one-h alo term ( Kashikawa et al. 200G: iLee et al 



20061 ). but is acceptable at low redshifts in the 



FUV since the one-halo term does not seem to 
play a major role as discussed sect. 13.21 

The mimimum masses of the DMH that pro- 
duce the bias derived for galaxies are lO 12 Af < 
M < 10 13 M Q at z > 2, 10 n M o < M < 1O 12 M 
at z ~ 1 and M < 1O 12 M at z < 0.4. There is an 
obvious degeneracy of the models at low redshifts, 
but the locally observed bias is definitely in the 
region of low cutoff masses. This is a hint that 
observed star-forming galaxies at low redshift re- 
side preferentially in less massive halos than high 
z star-forming galaxies. 



4.1. Redshift evolution of the bias 

Figure [TU] shows as symbols the redshift evo- 
lution of the bias parameter measured at 8 
h~ x Mpc defined as frs = vs,,g/&&,m (see e.g. 
M agl iocchetti et al. 2000) for the different sam- 



4.2. Bias and galaxy number density 

In fig. QT] we show the bias a function of the 
galaxy number density n ga i for UV-selected sam- 



4 For Adelberger et al. (2005) subsamples, 5 values are not 
available; we assumed that the slopes are the same than 
those of their global samples. The expected relative error 
on the inferred bias is lower than 10% if < S < 1. 
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Fig. 11. — Bias of restframe /TV-selected galaxies as 
a function of galaxy number density. The curves show 
the expected relation for the effective bias of Dark 
Matter halos at redshifts z = 0,0.5,1,2,3, from bot- 
tom to top. We plot only independent samples here. 
The legend is the same as in fig. [10] 



pies and the predicted relation between the effec- 
tive bias and the number density of DMHs at z = 
0, 0.5, 1, 2, 3 (curves from bottom to top). At high 
redshift (z > 1), we observe the well known lu- 
minosity segregation effect, brighter galaxies (less 
abundant) having a larger bias, in good agreement 
with DMH models predictions (the less abundant, 
the more clustered). In contrast at low z [z < 1), 
a significant departure to this relation is observed. 
At z ~ 0.9, the CFHTLS data show a bias slightly 
lower than the expected one according to the ob- 
served density with n ga i approximately 3 times 
lower than expected. This seems even worse for 
our brightest samples at z ~ 0.3, as these galaxies 
are about 10 times less numerous than expected 
according to their bias values. In the model dis- 
cussed here, we implicitly assume that one DMH 
hosts one galaxy, which provides a fairly reason- 
able description of the observations at high z, 
to the level of precision allowed here (see e.g. 
Adelberger et al.ll2005l ; lOuchi et alj|2004 for more 
detailed discussions on this point). At z < 1, this 
assumption may be not valid anymore and our re- 
sults suggest that star forming galaxies (especially 



the brightest) are not hosted by a significant frac- 
tion of the DMHs with similar clustering proper- 
ties. This implies that the DMH occupation frac- 
tion, that is roughly > 0.5 at high redshift (z > 2), 
drops to 0.3 and 0.1 at z = 1 and z — 0.3 respec- 
tively. 

4.3. Bias and FUV LD fraction 

The very limited overlap in FUV luminosities 
of the data at different redshifts does not allow 
a derivation of the bias evolution with redshift at 
fixed FUV luminosity. However, despite the fact 
that low z samples reach fainter luminosities than 
high z ones, they happen to probe the same frac- 
tion of total FUV luminosity densities, owing to 
the strong e volution of the FUV luminosity func- 
tion with z (jArnouts et alJl2005f ). In particular at 
all redshifts the samples are able to probe the bulk 
of star formation, i.e. they encompass a fraction 
of the FUV luminosity density (LD) greater than 
0.5. This can be seen in fig. [12] where we show 
the bias as a function of the fraction of the total 
FUV LD enclosed by the different samples. This 
favorable situation allows us to track the evolution 
with redshift of the clustering at a fixed fraction 
of the FUV LD, an essentially constant fraction of 
the star formation rate. 

The fraction of the LD for each sample is com- 
puted by comparing the total LD at the rele- 
vant redshift (from the FUV Lum i nosity func- 



tion parameters of lArnouts et alj (|2005i ) and 



ISawicki fe Thompsonl 12006)). to the LD enclosed 



by each sample according to its flux limits con- 
verted to FUV luminosity cutfl 
Note that we do not attempt to correct for galaxy 
internal dust attenuation. As brighter UV galax- 
i es present higher extinction in the local Universe 
( Buat et al]l2005h , the strong brighteni ng of FUV* 
with redshift (e.g. Arnouts et al. 2005! ) may intro- 
duce a small bias as a given LD fraction could 
not correspond exactly to the same star formation 
rate fraction at the different redshifts we explore 



5 Note that while LBG samples are by construction 
volume-limited, we cannot adopt this approach for the 
GALEX samples due to limited statistics. This means 
that faintest galaxies are underrepresented, especially in 
our higher redshift GALEX samples. However, as we do 
not observe a strong luminosity dependence of the bias 
within the GALEX samples, we expect this has only a 
small impact. 
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Fig. 12. — Bias of restframe (TV-selected galaxies 
(symbols as in figure II Op as a function of the frac- 
tion of the FUV total luminosity density represented 
by each sample. The legend is the same as in fig. [10] 



tocorrelation function well approximated by a 
power law in the range 0.2 — 5 Mpc. 

At z ~ 1, the correlation length of the rest- 
UV selected galaxies from CFHTLS u' data is 
found compa rable to those of the emission-line 
samples from ICoil et al.l (|2004) in the same red- 
sh ift range, but sl ightly higher than those obtained 
bv lMeneux et all {2006) for late type and irregular 
galaxies (their types 3 and 4, selected in the visi- 
ble). At those redshifts, according to our CFHTLS 
sample, star forming galaxies are modestly biased 
with (6g) = 1.26±0.06, which under the linear bias 
hypothesis implies they are closer tracers of the 
mass distribution than their higher redshift coun- 
terparts. As oppos ed to the dependence found 
at redshifts above 2 ( Giavalisco &: Dickinson|[200lt 



Foucaud et al.l 120031 ) , no strong positive correla- 
tion between the bias and the FUV luminosity is 
observed in the local universe, but rather a slight 
anti-correlation or no correlation. At z < 0.4, 
given the error bars, the mean bias is consistent 
with 0.8 for all GALEX samples ((b 8 ) = 0.79±£;^) 
independently of the UV luminosity. 



here. 

The plot confirms the result already apparent 
in fig. [TO] of a significant decrease of the bias of 
UV selected galaxies from high to low redshifts, 
but now selected on the basis of a physically de- 
fined parameter, the fraction of the FUV luminos- 
ity density. Near an LD fraction of 0.5, the bias 
is divided by a factor 3, between redshifts near 3, 
shifting from 2.5 in the redshift range 2 — 3 down 
to 0.8 in the local universe. 

5. Discussion 

In paper I we reported on the overall clus- 
tering properties of the UV-selected galaxies in 
GALEX samples, the largest ones available to 
date at low redshift and at these wavelengths. 
These samples allow for the first time an investi- 
gation of the clustering properties of UV-selected 
galaxies as a function of different parameters at 
z < 1, which can be compared to higher redshift 
samples also selected in the rest frame FUV. 

The measurements from the GALEX sam- 
ples confirm previous results for rest-UV selected 
galaxies at low re dshifts indicating that they are 
weakly clustered ( Heinis et al. 2004J), with an au- 



5.1. Migration of 'the bulk of star formation sites 
from z = 3 to the local universe 

In this study, we find a decrease by a factor 3.1 
of the bias with respect to mass, from redshifts 
near 3 to the local universe in the UV flux-limited 
samples, and more importantly in samples selected 
in UV luminosity so that they encompass a con- 
stant fraction of the luminosity density at all z. 
This decrease is sl ightly larger than th e factor 2.7 
derived from the iMo fc White] (|2002f ) model for 
the M > 10 12 Mq haloes that host most star for- 
mation at redshift 3, an indication that star form- 
ing galaxies tend to be hosted by haloes of lower 
mass in the local universe. This is the main con- 
clusion of the present study. 

The "downsizing" scenario (ICowie et al.lll996t 



Juneau et al.ll2005 ; Bundv et al.l2005f Heavens et al 



2004) states that the star formation shifts from 
high stellar mass systems at high redshift to low 
ones at low redshift. Our results extend this vi- 
sion in the sense that the same trend is observed 
for the mass of the dark matter halos that host 
actively star forming galaxies. 

The DMH mass migration of the bulk of 
the star formation might be associated with 
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regions of different densities. At high red- 
shifts, LBGs studies show that active star for- 
mation traced by the UV light resides preferen- 
tially in overdense regions (lAdelberger et al.l[l 998; 



Blai zot et all 12004 lGiavaliscol l2002l; ISteidel et al 
19981: iTasker fc Brvanl l2006h . At low redshift, 



Abbas fc Shethl ( 2005h showed that the slope of 
the fitted power law is steeper in underdense re- 
gions, and that the correlation length is smaller. 
The observed steeper ACFs for the more UV- 
luminous galaxies at low z suggest that the most 
star-forming objects reside preferentially in re- 
gions where the local galaxy density is lower than 
for the fainter ones, a result in agreement with 
direct optical based studies of star formation as 
a function of galaxy density in the l ocal universe 
([Gomez et al.ll2003HLewis et al.ll2002t) . 
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